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ABSTRACT 

The thermal stability of rotating, stratified, unmagnetized atmospheres is studied 
by means of linear-perturbation analysis, finding stability, overstability or instability, 
depending on the properties of the gas distribution, but also on the nature of the 
. perturbations. In the relevant case of distributions with outward-increasing specific 

| entropy and angular momentum, axisymmetric perturbations grow exponentially, un- 

less their wavelength is short enough that they are damped by thermal conduction; 
non-axisymmetric perturbations typically undergo overstable oscillations in the limit 
of zero conductivity, but are effectively stabilized by thermal conduction, provided 
rotation is differential. To the extent that the studied models are representative of the 
poorly constrained hot atmospheres of disc galaxies, these results imply that blob- 
' like, cool overdensities are unlikely to grow in galactic coronae, suggesting an external 

origin for the high-velocity clouds of the Milky Way. 

Key words: hydrodynamics, instabilities, ISM: kinematics and dynamics, ISM: 
clouds, galaxies: formation 
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1 1 INTRODUCTION 

CO ' Galaxy clusters and massive elliptical galaxies are embedded in hot atmospheres of virial-temperature gas revealed by X-ray 
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observations. Lower-mass galaxies are also believed to have hot gaseous coronae, difficult to detect because the gas is very 
rarefied. If the coronal gas were thermally unstable (Field 1965) it could fragment into cold gas clouds with substantial 
implications for the dynamics of cooling flows (Mathews & Bregman 1978; Cowie, Fabian, & Nulsen 1980; Nulsen 1986), 
but also for galaxy formation and for the origin of the Hi high- velocity clouds of the Milky Way (Mailer & Bullock 2004; 
Peek, Putman, & Sommer-Larsen 2008; Kaufmann et al. 2009). In fact, linear-perturbation analysis has shown that the X-ray 



S-h ' emitting hot atmospheres of massive elliptical galaxies and galaxy clusters are likely stabilized against thermal instability by 
9^ i a combination of buoyancy and thermal conduction (Malagoli, Rosner & Bodo 1987, hereafter MRB; Balbus & Soker 1989; 
Tribble 1989). The coronae of disc galaxies are expected to have lower gas temperatures and densities, but their detailed 
properties are poorly known, because — so far — they eluded detection in X-rays. Nevertheless, their physical parameters are 
not totally unconstrained: for some nearby massive disc galaxies upper limits to the total X-ray luminosity of the hot halos 
are available (Rasmussen et al. 2009), while in the special case of the Milky Way different pieces of information can be used 
to constrain the physical properties of the corona (e.g. Spitzer 1956; Fukugita & Peebles 2006). Studying model coronae 
consistent with these constraints, Binney, Nipoti & Fraternali (2009, hereafter BNF) showed that also the coronae of disc 
galaxies are likely thermally stabilized by buoyancy and thermal conduction, at least if rotation of the gas can be neglected. 
There are no significant constraints on the angular momentum of galactic coronae, but one might expect that they are 
characterized by at least slow rotation. It is well known that there is an important interplay between rotation and convection 
(e.g. Tassoul 1978; Balbus 2000, 2001), so rotation might influence the stabilizing effect of buoyancy and then the overall 
problem of the thermal stability of galactic coronae. 

This paper addresses the question of the thermal stability of rotating galactic coronae. For this purpose one needs to 
study the stability against non-axisymmetric perturbations of a differentially rotating, stratified gas in the presence of thermal 
conduction and radiative cooling. This is a special case of the general problem of the stability of a rotating stratified gas, which 
has been widely studied in the astrophysical literature, with specific applications to rotating stars (e.g. Cowling 1951; Goldreich 
& Schubert 1967; Lebovitz 1967; Fricke 1968; Lynden-Bell & Ostriker 1967; Tassoul 1978; Knobloch & Spruit 1982; Lifschitz 
& Lebovitz 1993; Menou, Balbus, & Spruit 2004) and accretion discs (e.g. Papaloizou & Pringle 1984; Balbus & Hawley 1991, 
1992; Ryu & Goodman 1992; Lin, Papaloizou, & Kley 1993; Brandenburg & Dintrans 2006). The stability-analysis techniques 
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developed in these previous studies can be applied to the problem addressed in the present paper. In particular, we will 
perform a linear analysis using Eulcrian perturbations, considering both axisymmetric and non-axisymmetric disturbances. 

It is reasonable to expect that galactic coronae are magnetized, though there are almost no constraints on the properties 
of their magnetic fields. For simplicity, in the present work we neglect magnetic fields, but we must bear in mind that they 
might have some influence on the thermal stability properties of the galactic hot gas, as suggested by previous studies of 
thermal instability in non-rotating cooling flows (Loewenstein 1990; Balbus 1991, see also BNF for a discussion). 

The paper is organized as follows. In Section 2 the problem is set up by recalling the relevant hydrodynamic equations. The 
thermal-stability analysis is then described for axisymmetric perturbations (Section 3) and for non-axisymmetric perturbations 
(Section 4). The implications for galactic coronae are discussed in Section 5; Section 6 concludes. A list of commonly used 
symbols is given for reference in Table 1. 



2 THERMAL INSTABILITY IN A ROTATING STRATIFIED GAS: GOVERNING EQUATIONS 



A rotating, stratified, unmagnetized atmosphere in the presence of cooling and thermal conduction is governed by the following 
equations for mass, momentum and energy conservation: 
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Here p, p, T and v are, respectively, the gas density, pressure, temperature and velocity, $ is the galactic gravitational 
potential, 7 = 5/3 is the ratio of principal specific heats, A is the cooling function, p is the mean gas particle mass in units 
of the proton mass m p and k is the thermal conductivity (Spitzer 1962). As we neglect magnetic fields, we treat the thermal 
conductivity as a scalar. In this approximation we can at most assume that the actual value of k is suppressed by tangled 
magnetic fields below Spitzer's benchmark value (e.g. Binney & Cowie 1981), but of course we cannot account for the effects 
related to the anisotropic heat transport expected in a magnetized medium (e.g. Balbus 2000; Quataert 2008). 

In cylindrical coordinates, assuming axisymmetric gravitational potential, the above hydrodynamic equations write: 
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The unperturbed atmosphere is assumed to be axisymmetric and close to hydrostatic and thermal equilibrium, in the sense 
that the system is approximately in a steady state over the time scales of interest even in the presence of radiative cooling and 
thermal conduction. Thus, the unperturbed fluid can be described by the time-independent axisymmetric pressure po, density 
po, temperature To and velocity field vo = (vor, vo z ,vo<f,), which satisfy equations (4-8) with vanishing partial derivatives with 
respect to t. The velocity components vqr and vo z are allowed to be non-null, because an inflow of gas can occur if cooling is 
not perfectly balanced by thermal conduction. Such a steady-state configuration is expected to be a reasonable approximation 
for the corona sufficiently far from the disc, where we are interested in studying the thermal instability. The model is not 
intended to apply near the disc and close to the galactic centre, where the interactions of the coronal gas with the star-forming 
disc and the central super-massive black hole are important. 

In the most general case the atmosphere is assumed to be differentially rotating, with angular velocity fl = votp/R 
depending on both R and z, but in the analysis it is convenient to distinguish cases with d£l/dz = from cases with 
dQ/dz / 0. The Poincare-Wavre theorem (e.g. Tassoul 1978) states that the surfaces of constant pressure and constant 
density coincide if and only if dQ/dz — 0: as a consequence, distributions with dQ/dz — are said barotropic (pressure is a 
function of only density), while distributions with dfl/dz / are said baroclimc (pressure is not a function of only density). 



Thermal instability in rotating galactic coronae 3 



Table 1. List of commonly used symbols. 

A p r = (dpo/9R)/p 
A pz = (dp /dz)/p 
A pR = (dp /dK)/p 
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Inverse of the pressure scale-length 
Inverse of the pressure scale-height 
Inverse of the density scale-length 
Inverse of the density scale-height 

Isothermal sound speed, normalized isothermal sound speed 
Derivative along surfaces of constant wave-phase 
Perturbation wave-vector 
Perturbation wave-number 

Normalized wave-vector components and wave-number 
Dopplcr-shiftcd perturbation frequency 

Pressure perturbation throughout the paper (pressure in Section 2) 
Unperturbed pressure 
Unperturbed specific entropy 

Temperature perturbation throughout the paper (temperature in Section 2) 
Unperturbed temperature 
Unperturbed velocity 

Components of the velocity perturbation throughout the paper (of the velocity in Section 2) 

Normalized components of the velocity perturbation 

Ratio of principal specific heats 

Local polytropic index of barotropic distributions. 

Local R and z logarithmic slopes of the unperturbed pressure 

Local R and z logarithmic slopes of the unperturbed density 

Local R and z logarithmic slopes of the angular velocity 

Thermal conductivity 

Radiative cooling function 

Mean gas particle mass in units of the proton mass m p 

Density perturbation throughout the paper (density in Section 2) 

Unperturbed density 

Normalized density perturbation 

Normalized time 

Perturbation frequency 

Dopplcr-shiftcd perturbation frequency 

Brunt- Vaisala, (or buoyancy) frequency (equations 18 and 21) 
Characteristic frequency associated with rotation (equations 17 and 33) 
Thermal-instability frequency (equation 15) 
Thermal-conduction frequency (equation 14) 
Characteristic frequency of dissipative processes 
Angular velocity 

Local R and 2 derivatives of f o<£ = QR 



To address the question of the thermal stability of the rotating corona, we study the behaviour of the system in the 
presence of small (linear) thermal perturbations. These perturbations are not expected to have particular symmetries in a 
real system, so the relevant case is that of non-axisymmetric perturbations in the presence of differential rotation. It is well 
known that in a differentially rotating system such a problem is complicated by the effect of the shearing background on the 
perturbations (Cowling 1951; Goldreich & Lynden-Bell 1965; Lynden-Bell & Ostriker 1967). The analysis is much simpler in 
the case of axisymmetric perturbations or uniform rotation: though these assumptions are not expected to apply in general 
to galactic coronae, exploring these cases is useful to understand the behaviour of the more realistic, but far more complex 
case. Therefore, before addressing the full problem of non-axisymmetric perturbations in a differentially rotating corona, 
we will consider the cases of axisymmetric perturbations in a differentially rotating atmosphere and of non-axisymmetric 
perturbations under the assumption of uniform rotation. 



3 THERMAL-STABILITY ANALYSIS: AXISYMMETRIC PERTURBATIONS 
3.1 Derivation of the dispersion relation 

Here the fluid is assumed to rotate differentially with O = z) and that the perturbations are axisymmetric. Linearizing 
the hydrodynamic equations (4-8) with perturbations of the form Fo + F exp(— iojt + ik R R + ik z z) (where Fo is the unperturbed 
quantity and |F| <C |-Fo|), in the limit of short- wavelength, low-frequency perturbations we get 



-iuip + ik R v R po + ik z v z po = 0, 
-iQv R p = -ik R p + A pR clp + 2f2«0po, 
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where o) = aj — k-vo=w — (kRVQR + k z vo z ) is the Doppler-shifted frequency, C 2 , = Po/po is the isothermal sound speed squared, 
= (dpo/dR)/po and ^4 pz = (dpo/dz)/po are the inverse of the pressure scale-length and scale-height, respectively. The 
following frequencies have also been defined: Qr = d(flR)/dR, Q, z = d(£lR)/dz, the thermal-conduction frequency 
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and the thermal-instability frequency 
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In terms of the defined quantities, the assumption of short-wavelength perturbations gives \k z \ 3> \ A pz \, \ j 

and Q, 2 , Q, R , Q, 2 <§; Cq/c 2 , while the assumption of low- frequency perturbations gives J 2 -C c\k 2 . In deriving the energy equa- 
tion (13) we used poT ~ —Top, because it can be shown that in the considered limit p/po <C p/po (see Appendix A). The 
system of equations (9-13) can be reduced to the following dispersion relation 1 for n = —icb: 



n + n uj d + (uj bv + Lu Iot )n + w Iot uj d = 0, 

where u} d = u} th + lu c is the characteristic frequency of dissipative processes, 
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InpoPo 7 is the unperturbed specific entropy and, following Balbus (1995), we introduced the 



(19) 



is the buoyancy term, so 
differential operator 

JL 
dR' 

which can be seen as taking derivatives along surfaces of constant wave phase. It may be useful to note that, by definition, 
w d > 0, while w| v and uj 2 ot are not necessarily non-negative, and that n is defined so that the perturbation evolves in time 
as oc exp(ni), thus the stable modes are those with Re(n) < 0. 



k z dz 



3.2 Limiting cases 

Before analyzing the dispersion relation (16) in the general case, it is useful to discuss some limiting cases, which can be 
obtained when one or two among the three characteristic frequencies ujbv , <^rot and oj d are zero. Let us start from the simplest 
cases, in which only one of these is non-null. 

(i) Case with Wd — and u I0 t = 0. In this case there is no dissipation and the fluid is either non-rotating (fi = 0) or 
rotating differentially with vanishing gradient of the specific angular momentum [d(Q,R 2 )/dR — when Q, = Q(R)]. From 
equation (16) we have the dispersion relation 

n = -wlv, (20) 

so we have stability if the square of the Brunt- Vaisala frequency w| v > 0. When Q, = Q(R) the Brunt- Vaisala frequency 
squared con be written 

.2 _ k 2 z c%A 2 vz ( 7 A fk R A p r 



•^-^(^-VKsS) • (21) 

where, using the fact that po is a function of only p (because the distribution is barotropic; see Section 2), we defined 



1 The same dispersion relation is obtained by working from the beginning in the Boussinesq approximation, i.e. neglecting the term 
— icjp in equation (9) and the term —iuip/po in equation (13). 
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V-Jr^. (22) 

din p 

which can be considered a local polytropic index. It follows that, independently of the value of kn/k z , the condition for 
convective stability is 7' < 7, i.e. Schwarzschild's criterion (e.g. Tassoul 1978). 

(ii) Case with cj B v = and cj ro t = 0. Here we assume that the characteristic frequencies associated with rotation and 
buoyancy are null. In the barotropic case these conditions are met when 7' = 7 (i.e. the radial gradient of the specific entropy 
is zero) and d(QR 2 )/dR — (i.e. the radial gradient of the specific angular momentum is zero). The dispersion relation is 

n — — LJd, (23) 

so we have thermal instability if ujd < 0, which is just Field's instability criterion (Field 1965). From the definition of 
cj d = Wth + w c , it is clear that the condition for thermal instability is that the growth rate of the thermal perturbation 
(|wth|) must be faster than conductive damping (we recall that uj c > 0, while typically wth < 0). For fixed unperturbed 
gas temperature To and pressure po, oj c increases for increasing perturbation wave- number k (equation 14), while uj t h is 
independent of k (equation 15), so there is a critical perturbation wavelength such that ujd < for longer wavelengths and 
u>d > for shorter wavelengths. 

(iii) Case with wbv = and cjd = 0. In the absence of buoyancy and dissipation we obtain the dispersion relation 

n = -Wrot, (24) 

so the stability criterion is w r 2 ot > 0. The value of u) 2 ot depends on the ratio kn/k z : w r 2 ot is positive for all values of kn/k z if and 
only if dft/dz = and d(Q,R 2 )/dR > 0, i.e. the specific angular momentum must increase outwards: this is just Rayleigh's 
stability criterion (Chandrasekhar 1961). 

In the following three limiting cases of the dispersion relation (16) only one among ujbv, <^rot and oj d is null. 

(i) Case with u; B v = 0. In the absence of buoyancy, but for uj TOt 7^ and uj d 7^ 0, the dispersion relation is 

(n 2 +w r 2 ot )(w d + n) =0, (25) 

which is just a combination of the criteria obtained in the points ii) and iii) above, so (when ujbv = 0) the presence of a 
gradient of the specific angular momentum does not modify the thermal-instability criterion in an interesting way. Specifically, 
when cjd < the medium is thermally unstable, independent of the presence and properties of rotation, while rotation can 
destabilize an otherwise thermally stable medium (wa > 0) if uj 2 ot < 0. Thus, the condition for stability is cod > and u) 2 ot > 0, 
which holds for all values of kn/k z if and only if dQ,/dz = and d(flR 2 ) / dR > 0. 

(ii) Case with cod = 0. In the absence of dissipation, but for w ro t =fc and wbv =fc 0, one obtains the dispersion relation 

n = -(«iv +Wrot). (26) 

which leads to the well-known convective stability criterion for a rotating, stratified fluid uj^+lu 2 ^ > 0, showing the stabilizing 
effect of rotation against convection. The inequality oj| v + Ld 2 ot > is verified for all values of ku/k z if and only if 

1 „ 1 dR 4 Q 2 „ 

Vp .Vso + ^3-^- >0 27 

7Po R* oR 

and 

d Po ( dR 4 n 2 ds dR 4 n 2 ds * n 



dz \ dR dz dz OR 

which is the classical Solberg-H0iland criterion (see Solberg 1936; H0iland 1941; Goldreich & Schubert 1967; Tassoul 1978; 
Balbus 1995). 

(iii) Case with ui I0 t = 0. When the fluid does not rotate, or, more generally, has a vanishing gradient of the specific angular 
momentum [d(QR 2 )/dR = in the barotropic case], but u>d ^ and uibv ^ 0, the dispersion relation is 

n 2 + uj d ri + cj| v = 0. (29) 

If the system is spherically symmetric, equation (29) reduces to the dispersion relation derived in MRB (see also BNF). When 
> we have stability (damping by thermal conduction) if w| v > 0, while we have convective instability if w| v < 0. When 
uJd < we have overstability if uj BV > |, thermal instability if < w BV < \ and convective instability if w BV < 0. 

3.3 Analysis of the general form of the dispersion relation 

Here we consider the dispersion relation (16) in the general case in which all coefficients are non-null. Dividing by |cjd| 3 we 
get 
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h 3 +h 2 + (w| v + £b 2 ot )h + cD r 2 ot = 0, if cu d > 0, (30) 
and 

n A - n + (w|v + Wr 0t )fi - cD 2 ot = 0, if w d < 0, (31) 

where the dimensionless quantities h = n/|wd|, w| v = u)% v /lo\ and <I) 2 ot = u 2 ot /u^ have been introduced. In both cases the 
discriminant of the cubic equation is 

A = -27a: 2 + (36y + 4)x - 4y(y 2 + 2y + l), (32) 

where x = w| v and y = w| v + uj 2 ot . 

Let us discuss first the case cjd > (equation 30 and diagram in the left-hand panel of Fig. 1). Applying the Routh- 
Hurwitz theorem (see Appendix B) we have that the real parts of all roots are negative (stability) if and only if w| v > 
and cJ 2 ot > 0. Thus the first quadrant of the diagram in the left-hand panel of Fig. 1 is a locus of stable configurations. 
Configurations in the other three quadrants can be either unstable or overstable. The blue curves in the diagrams correspond 
to discriminant A = 0. In the bottom-left area defined by these curves A > 0, thus we have three real roots (at least one 
negative), so the corresponding configurations are unstable. In the other regions A < 0, so we have one real root (hi) and two 
complex conjugate roots (fi2 and h$). The sign of the real root hi can be determined because we know that hihihz = — u) 2 ot 
(see Appendix B), and hzhs is obviously positive. Thus we have instability in the upper-left region and overstability in the 
bottom-right region. In summary, the condition for stability against axisymmetric perturbations is wjlv > and cj 2 ot > 0, 
which holds for all values of fcjj/fc z if dfl/dz = [thus Q — Q,(R) and po — po(po)], d(QR 2 )/dR > and dlnpo/dlnpo < 7- 
We note that the analysis of the case Wd > was carried out also by Lifschitz & Lebovitz (1993) in the context of the study 
of the stability of rotating stars in the presence of radiative diffusion (see also Goldreich & Schubert 1967; Sung 1974b, 1975; 
Balbus 2001). 

Let us move to the case oj d < (equation 31 and diagram in the right-hand panel of Fig. 1). The same line of reasoning 
as for the case uj d > leads to the following conclusions. Given that the coefficient of h 2 is negative, we know from the 
Routh-Hurwitz theorem that in no case all the real parts of the roots are negative. In other words, all configurations will be 
either unstable or overstable. The blue curves in the diagram in the right-hand panel of Fig. 1 correspond to discriminant 
A = 0. In the bottom-left area defined by these curves A > 0, thus we have three real roots (at least one negative), so the 
corresponding configurations are unstable. In the other regions A < 0, so we have one real root and two complex conjugate 
roots: in this case the sign of the real root is the same as the sign of <I> 2 ot . Thus we have overstability in the upper-left region 
and instability elsewhere. It follows that a necessary condition for overstability is w| v > and uj 2 ot < 0, which holds for all 
values of kn/k z if dQ/dz = [thus = Q(R) and po = po(po)], d(QR 2 )/dR < and 7' < 7. Note that this condition is not 
sufficient, because for positive, but small enough values of cjj|v instability replaces overstability even if cj 2 ot < (see top-left 
quadrant of the right-hand panel of Fig. 1). The non-rotating case (briefly discussed at the end Section 3.2) is obtained from 
the diagram at uj 2 ot = (overstability if ti| v > |, and instability at lower values of w| v ; see MRB). 



3.4 Discussion 

The above results indicate that the signs of cjj|v an d w ?ot ar e crucial for the thermal-stability properties of a configuration. 
Both wlv an d w ?ot depend not only on the properties of the unperturbed distribution, but also on the wave-vector of the 
perturbations through the coefficient fc 2 /fc 2 , and through kn/k z , which appears in the differential operator T> defined in 
equation (19). A general result, valid for both positive and negative cj d , is that if fi depends on z the system is unstable, in 
the sense that it is always possible to choose kn/k z such that we are in the instability region of parameters. To look for stable 
or overstable configurations we must consider barotropic distributions, characterized by = Q,(R) and po = po(po)- In this 
case uj 2 ot is given by 

A (33) 
so it reduces to (k z /k) 2 times Rayleigh's discriminant (Chandrasekhar 1961). The condition for co 2 ot > is 
d In Q 

din* > ^ ( 34 ) 
which is Rayleigh's stability criterion, requiring that the specific angular momentum must increase outwards. In addition, for 
barotropes the Brunt- Vaisala, frequency is given by equation (21), so the condition for Wg V > is Schwarzschild's criterion 

7 < 7- (35) 

In most applications we are interested in the case of convectively stable systems (satisfying condition 35) with outward- 
increasing specific angular momentum (satisfying condition 34). Coming back to the diagrams in Fig. 1 we conclude that 
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Figure 1. Domain of stability (red), overstability (green) and instability (blue) against axisymmetric perturbations in the plane iigy = 
^Bv/i^dl 2 versus k'rot = aJ ?ot/l a; d| 2 i f°r a differentially rotating, stratified fluid, when > (left-hand panel) and u>^ < (right-hand 
panel). uj^ ot and <^g V are given, respectively by equations (17) and (18). 




in the presence of cooling and thermal conduction, and assuming axisymmetric perturbations, these systems are stable if 
ojd > (damping by thermal conduction), while are unstable if tv^ < 0, independent of the value of the ratio lo% v /ui^ and 
on of the value of to^ ot > 0. This means that in the case a>d < (i.e., when the perturbation wavelength is short enough that 
the thermal-instability frequency is higher than the conductive-damping frequency) even arbitrarily slow rotation changes 
qualitatively the behaviour of the system transforming overstability (found in the absence of rotation; MRB) into instability. 
Physically, this different behaviour is due to the fact that, when the perturbations are axisymmetric, conservation of angular 
momentum tends to inhibit convection (Cowling 1951) and thus to prevent overdense regions from falling inward far enough to 
oscillate around the radius proper to their particular value of specific entropy. Note that the transition between overstability 
and instability is not abrupt, because the growth rate of the unstable mode goes to zero when w^ot ~* (for oi BV > j). This 
is apparent from Fig. 2, plotting the instability growth rate n, which is the real root obtained by solving equation (31) for 
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cj, = cj,,+cj >0 cj. = cj., +cy <0 

d tn c d Lh c 




10 20 30 40 10 20 30 40 
"vol S ?ot 

Figure 3. Domain of stability (red), overstability (green) and instability (blue) against non-axisymmetric perturbations in the plane 
(I)g V = cjgy/loidl 2 (equation 45) versus ui^ ot = ui^ ot /\u^\ 2 (equation 46), for a uniformly rotating, stratified gas, when cjj > (left-hand 
panel) and when w d < (right-hand panel). The arrows show the effect of increasing the azimuthal wave- number k^ = m/R from = 
(filled circles) to \k^/k z \ 3> 1, for representative cases with 7' = 1 (solid lines) and 7' = 3 (long dashed line), assuming k z /kn = 2, 
|u> d |/fi = 0.3, c = 2, T pz = -1.5, and F pR = -1.25. 

^rot > and [Dbv > 0, as a function of uj^ ot and cDbv- The growth rate, at fixed w B v, increases gradually from for ui^ ot 
increasing from 0; the larger uj bv the larger must be uj^ ot to give substantial growth rates. 



4 THERMAL-STABILITY ANALYSIS: NON-AXISYMMETRIC PERTURBATIONS 

In the present Section we move to the next step of the stability analysis, which is the study of non-axisymmetric perturbations. 
Formally, it would be sufficient to perform the non-axisymmetric perturbation analysis for those configurations that are stable 
or overstable against non-axisymmetric perturbations, which might turn out to be unstable if more general perturbations are 
considered. Nevertheless, we will perform the non-axisymmetric perturbation analysis also for configurations found unstable 
against non-axisymmetric perturbations, because it is interesting to verify whether axisymmetry is a necessary condition for 
a perturbation to grow. 

There are reasons to expect that non-axisymmetric thermal perturbations behave differently from axisymmetric ones. 
Cowling (1951) showed that rotation has a stabilizing effect against convection for axisymmetric perturbations (a consequence 
of angular momentum conservation), but not necessarily for non-axisymmetric perturbations: in the latter case the stabilizing 
effect of rotation is less and it vanishes for specific perturbations. This suggests that rotating galactic coronae (in the absence 
of heat conduction) might be thermally unstable against axisymmetric perturbations, but not against non-axisymmetric 
perturbations, because buoyancy should prevail, provided the gas entropy increases outwards (see also Nipoti 2010). 

We now address quantitatively the problem of stability against non-axisymmetric perturbations. As mentioned above, in 
this case the analysis is rather complicated if the corona rotates differentially (Goldreich & Lynden-Bell 1965). Therefore, we 
first try to get an estimate of the effect of the deviation from axisymmetry of the perturbations by considering the simpler case 
of non-axisymmetric perturbations in a uniformly rotating corona (Section 4.1), and we defer the treatment of the differentially 
rotating case to Section 4.2. 

4.1 Uniform rotation 

Linearizing the hydrodynamic equations (4-8) with perturbations of the form Fo + Fexp(— iuit + ikaR + ik z z + im(j>), where 
Fo is the unperturbed quantity and \F\ <t. \Fo\, and assuming angular velocity SI independent of position, we get: 

— \ujp + iknVRpo + ik z v z po + ik^v^po = 0, (36) 
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-iuvRpo = -\k R p + A pR c p + 2Qv. 
-iQv z po = -ik z p + A pz c 2 -,p, 
-iCjv^po + 2v R p Q, = -ifc^p, 



P P 

iw h ijcj h v R (A pR - yApR) + v z (A pz - jA p 

Po Po 



(37) 
(38) 
(39) 

(40) 



where now k 2 = k R + k z +k 2 and Co = cj — k-vo = uj — (kRVoR + k z v 0z + k^v ^) , with k<p = m/R. In deriving the equations above 
we have assumed, as done in Section 3, short wavelengths {\k^\, \kii\, \k z \ >• |-A p .r|, \A pz \, \A p r\, \A pz \, fi/co), low frequencies 
(cu 2 -C Cofc 2 ) and poT = —Top. The system of equations can be reduced to the dispersion relation 



n A + n 2 io d + (w|v + w 2 ot )n + io 2 ot co d 



0. 



where now 



2 



7Po k 2 



-n -n i fdPodso dp ds \ 



; 2 Ji2 
fc 2 7 



7' 1 



+ _T I 1 + 



(41) 



(42) 



is the buoyancy term, which for = reduces to the Brunt- Vaisala, frequency in equation (21), with 7' given by equation (22), 
and 



2 _ k z 2 
w rot = ^4fi 



(43) 



(44) 



is the rotation term, which is just io 2 ot as defined in equation (17) for position-independent fi. Note that in deriving the 
dispersion relation (41) we exploited the relation 

<9po <9so _ dpo ds 
~dz~~dR ~ ~dR~dz~' 

which holds because the surfaces of constant pressure and constant specific entropy coincide in fluids rotating with f2 inde- 
pendent of z (i.e. with barotropic distributions; see Section 2). 

The dispersion relation (41) for non-axisymmetric perturbations in a uniformly rotating medium is formally identical to 
the dispersion relation (16) for axisymmetric perturbations. We just note that now ajj| v depends on (equation 42) and 



> (equation 43). When ui d = we obtain the dispersion relation n 



(ui BV + w rot ) (with lj bv and to rot given by 



equations 42 and 43), so the necessary and sufficient condition for convective stability against non-axisymmetric perturbations 
is oj bv + io Iot > 0, showing the stabilizing effect of uniform rotation against convection (in accordance with previous studies; 
Cowling 1951; Sung 1974a, 1975; Ryu & Goodman 1992). When LOd 7^ 0, the analysis of the dispersion relation (41) is the 
same as that of Section 3 and leads to the diagrams reported in Fig. 3, where now 



2 

-2 _ ^BV 
W BV = — — 
W d 



k 2 z c r pz ( 7 



and 



2 

^rot 



k 2 7(1)2 V 7' 



J_kj 

k 2 '' 



k_R 

k z 



r 



k 2 
k 2 



1 + 



1 pR 

r 2 
j- pz 



^d ^d 

having introduced the dimensionless quantities Co 



co/QR, uJd = Wd/fi, 



pR 



d In po 
dlaR 



= RA 



pR 



and 



_R dhip 
lpz ~ z 9 In Id - KA ^- 



(45) 



(46) 



(47) 



For both uid > (left-hand panel in Fig. 3) and ojd < (right-hand panel in Fig. 3) the domain of stability, instability and 
overstability are the same as in the first and fourth quadrants of the corresponding diagrams of Fig. 1. 

The deviation from axisymmetry of the perturbation is measured by the ratio k^/k 2 , on which both uj 2 ot and cj| v depend. 
When the properties of the unperturbed system (in practice, Cod, Co, Y p r, T pz and 7') and the ratio of the R and z wave- 
numbers kR/k z are fixed, the position in the Cj 2 ot -w BV plane depends only on k^/k 2 . The sign of cj| v depends only on 7' 
(equation 42), so the quadrant in which a system is located does not change by varying the ratio k^/k 2 at fixed 7'. Therefore, 
when iOd > a uniformly rotating system that is thermally stable against axisymmetric perturbations is stable also against 
non-axisymmetric perturbations (see also Sung 1975). Similarly, when Lo d < a uniformly rotating system that is thermally 
unstable against axisymmetric perturbations is unstable also against non-axisymmetric perturbations. 

In order to illustrate the effect of varying k^/k 2 , in Fig. 3 we plot the results for representative cases with 7' < 7 (namely 
7' = 1; solid arrows) and 7' > 7 (namely 7' = 3; long dashed arrows), assuming k z /kR = 2, |a>d|/fi = 0.3, Co = 2, Y pz = —1.5, 
and YpR — —1.25. Of course the specific behaviour of the system depends on the choice of the parameters, however the 
examples shown in Fig. 3 represent qualitatively the trends of typical models in the four quadrants. When all the other 



10 C. Nipoti 



parameters are fixed, the modulus of ajj| v increases for increasing k 2 /k 2 . Therefore, when the system is convectively stable 
(w| v > 0) increasing k$ (i.e., considering non-axisymmetric perturbations with larger m or shorter azimuthal wavelength) 
moves our system towards the top-left of the first quadrants of the diagrams in Fig. 3. It is interesting to note that when 
\ktf,/k z \ — > oo, the system can become overstable, because uj 2 ot — > (we recall that for aj<j < the system is overstable if 
Wrot = and ti| v > \). Systems with w| v < move, for increasing k<f,, towards the bottom-left part of the diagrams in Fig. 3: 
in these cases, non-axisymmetric effects cannot change the unstable nature of a system with cjd < 0, but it is interesting to note 
that a system with ojd > 0, overstable against axisymmetric perturbations, can be unstable against high-m non-axisymmetric 
modes. 



4.2 Differential rotation 

Here we consider the case of non-axisymmetric perturbations in a differentially rotating corona, which is complicated by the 
effect of the shear on the perturbations. A possibility would be to consider disturbances in the form Fq + F(R, z) exp(— iuit + 
im<f>), where Fo is the unperturbed quantity and \F\ <C \Fo\ (see e.g. Cowling 1951; Lin & Shu 1964; Sung 1974a, 1975; Fujimoto 
1987; Hanawa 1987; Bertin et al. 1989): linearizing the hydrodynamic equations with these perturbations leads in general to 
a system of partial differential equations. As an alternative, one can adopt shearing coordinates (Goldreich & Lynden-Bell 
1965) to reduce the problem to a system of ordinary differential equations (see also Balbus & Hawley 1992; Brandenburg & 
Dintrans 2006; Balbus et al. 2009). Here we follow the latter approach, and in particular we adopt the formalism used by 
Balbus & Hawley (1992) in their study of non-axisymmetric perturbations in a differentially rotating magnetized disc. 

Let us consider an unperturbed distribution that is a solution of equations (4-8) with vanishing time partial derivatives 
and vo = (0, 0, vo</>), where V04, = fl(R, z)R (for simplicity we assume here vor = vo z = 0). Following Goldreich & Lynden-Bell 
(1965) we perform the change of coordinates <f>' — (f> — Q(R,z)t, R' = R, z — z and t' = t. The partial derivatives are 
transformed as follows: 

I = w m 

OR OR' dRdf 1 ; 

d__d__ t dV_d_ 

dz dz> dz d<t>' ' { ' 

3D 3 

i* = 8F (51) 

We perturb the system with small non-axisymmetric disturbances. In the primed coordinates the perturbation takes the form 
of a plane wave, so we can write the perturbed quantities as Fo + F(t') exp(ik R R' + \k' z z + imcj)), with k' R , k' z and m constant. 
It is convenient to work in the Boussinesq approximation: linearizing equations (4-8) we get 

k R {t')v R + k z (t')v z + fe^v^ = 0, (52) 
-^j-Po + ik R (t')p - ApRCgp - 2nv 4 ,p = 0, (53) 

^fpo + ik z (t')p- A pz clp = 0, (54) 

^rn-po + v R p (Cl + Q R ) + v z p Q, z + ik^p = 0, (55) 
at 

-1^£:— + vr(A pR -iA pR )+v z (A pz -iA pz ) -7— Mt')+ w th] = 0, (56) 
at' pa pa 

where k^ = m/R' , 

kR(t') = k'R-mt'^ (57) 
and 

kJt') = k' z -mt'^. (58) 
oz 

Note that now cj c , still defined as in equation (14), depends on time, because it depends on k 2 (t') = fcij(t') 2 + k z (t') 2 + k\. 
As implementations of the Boussinesq approximation we neglected the term proportional to dp/dt' in equation (52) and the 
term proportional to dp/dt' in equation (56), and we assumed poT — —Top. 

The system of five ordinary differential equations (52-56) in the unknown vr, v z , V4,, p and p can be simplified by 
eliminating v$ and p. In particular, in the momentum equations, v<p can be eliminated by using the mass conservation 
equation (and its time derivative): 



Thermal instability in rotating galactic coronae 11 



k R (t') k z (t') 

cky _ k R (t') dv R k z (t') dv z . . 

dt' ~ fc dt' fe^, dt' ' 1 j 

The pressure perturbation p can be eliminated by combining the momentum equations as follows: we subtract the z equation 
(multiplied by k R ) from the R equation (multiplied by k z ), and we subtract the z equation (multiplied by k^) from the <f> 
equation (multiplied by k z ). After rearrangement and simplification, we end up with the following system of three coupled 
ordinary differential equations (in the variables vr, v z and p), which fully describes the evolution of the non-axisymmetric 
perturbations in a differentially rotating corona: 

dv R 2k R k 4> Q, / 9 lnQ k 2 z \ 2k R k 4 ,n ( R dlxiQ, k z fc|+_fc| \ c\A vR I ' k\ + k\ A pz k R k z \ 

dt ~ k 2 \dlnR kl) VR+ k 2 \zdln\z\ k R k 2 J Vz + p„ ^ k 2 A pR k 2 J P ' ( ' 

dv z _ 2fc z fc n ( d ln„ k 2 R + kl \ '2k z k^l ( R dlnQ. k z k R \ , (? A pz ( k\ + k% A pR k R k z \ 

dt ~ k 2 \dlnR + k% ) VR+ k 2 \zd\n\z\ + k 2 ) Vz + p { k 2 A pz k 2 J P ' 

^ = — (A P R - "/A p r) vr + — (A pz — "fA pz ) v z - („ c + „th) p, (63) 
dt 7 ' 7 

where we dropped the explicit time dependence of k R , k z and uj c . It is useful to rewrite the system above more compactly 
using only dimensionless quantities: 



dv R 2k R , r2 \ _ , 2k R 

dr k 2 K ' k 2 



r S2z -i_-(fc 2 + l) 



2,- I k z + 1 V pz kRk z 1 

^ +Co r pfl (-__-_-_-]p, (64) 



d7 = 17 ( rM + fc « + : ) VR + 17 ( r ^ + KkR ) v * + c « r ^ (^-p- -rv^2-J p ' (65) 

^- = — (TpR - ^T p r)v r + — (F pz - "/T pz )v z - (cD c fc 2 + _th)p, (66) 

QT 7 7 

where % = vr/QR, v z = v z /QR, p = p/po, r = t'Q, k R = k R /k^, k z = k z /k 4> , k = k/k^, 

_ dlnpp _ _ R dlnpo _ 

rpR = ^R- RApR > rpz = ~dh[\ Z ~\- RApz > (67) 

_ dinn _ i? ainn 

a In Jt z a m 1 2 

r p _R and r pz are defined in equations (47), 

. _ u c kl i / 7 -i\ fc^r 7/2 

np =n( v -rj^— ' (69) 

and 

-th = (70) 

In equations (64-66), besides the unknown vr(t), v z (t) and p(r), the only other time-dependent quantities are 

k R (r) = k' R - T sm r, (71) 

k z (r) = k' z -T nz T, (72) 

and fc 2 (r) = k\ + k z + 1, where fc^. = k' R /k^ and fc z = k' z /k^ (note that <D C does not depend on time). The system of 
ordinary differential equations is thus completed by specifying at r = the values of vr, v z and p, and by choosing the 
values of the following parameters: k' R , k' z , Fqr, Tn z , T p r, Y pz , Y p r, T pz , Co, Co c and „ t h- A full exploration of the parameter 
space is prohibitive, so we will present solutions for specific relevant cases, obtained by numerically integrating the system of 
equations (64-66) with a fourth-order Runge-Kutta method. 



4-2.1 Barotropic distributions 

Let us focus first on solutions of the system of equations (64-66) in the case of barotropic distributions (see Section 2), for 
which Q — n(7?) and po = po(po), with local polytropic index 7' = dlnpo/dlnpo. We are interested in estimating the effect 
of the non-axisymmetry of the perturbations and of the differential rotation on the thermal stability of the fluid, so it is 
convenient to make a comparison with the results obtained for axisymmetric perturbations (Section 3) and uniform rotation 
(Section 4.1). In those cases, we found it useful to explore separately the cases tud = 0, u)d > and cjd < 0, where u)d = uj t h + ^c 
is the characteristic frequency of dissipative processes. In the case of non-axisymmetric perturbations and differential rotation 
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cj d = 0,{Co c k 2 +cD t h) depends on time through fc(t'), so for comparison with the previous analysis it is convenient to distinguish 
the following cases: 

(i) Case with Hi c — and w t h = 0. This is the non-dissipative case, which has been studied in previous work on convective 
instability in rotating stratified fluids (Cowling 1951 ; Sung 1975). In accordance with Sung (1975), we find that the combination 
of differential rotation and non-axisymmetry of the perturbations can have a destabilizing effect against convection. For 
instance, we find that — for sufficiently large values of kj, — adiabatic (7' = 7) configurations are stable for Tqr = (uniform 
rotation), but are unstable for Fqr / (differential rotation), even when the specific angular momentum increases outwards 
(T (m > -2). 

(ii) Case with cD c > and w t h = 0. This is a dissipative case, with positive dissipative term, analogous to the model 
with radiative diffusion studied by Sung (1975). For axisymmetric perturbations or uniform rotation, when ojd > the only 
stable configurations are those characterized by 7' < 7 and d(QR 2 )/dR > (i.e. T^r > —2); the other configurations are 
unstable or overstable (see left-hand panels of Figs. 1 and 3). We calculated numerical solutions for configurations with 7' < 7 
and Tqe > —2 (known to be stable against axisymmetric perturbations), finding them stable also against non-axisymmetric 
perturbations for a wide range of values of k,R and k z (see also Sung 1975). 

(iii) Case with cD c = and uj th < 0. This is a dissipative case, with negative dissipative term, to be compared with the 
cases of axisymmetric perturbations or uniform rotation with a>d < (right-hand panels of Figs. 1 and 3). In those cases, the 
systems were found unstable in the relevant parameter regime 7' < 7 and TnR > —2 (outward-increasing specific entropy 
and angular momentum). We explored the stability against non-axisymmetric perturbations for several configurations with 
7' < 7 and F(ir > —2 finding that the combination of non-axisymmetry of the perturbations and differential rotation can 
turn instability into overstability. This can be seen from Fig. 4, showing, for example, the time evolution of non-axisymmetric 
perturbations in three configurations (models 4a, 4b and 4c) differing only for the values of Fur (measuring the degree of 
differential rotation), k' R and k' z (measuring the deviation from axisymmetry of the perturbations). The modes considered in 
models 4a and 4b have k^ of the order of kR and k z , so they are relevant to blob- like perturbations. In model 4a rotation 
is uniform, and the system is unstable, consistent with the calculations of Section 4.1. Model 4b differs from model 4a only 
because differential rotation (with Tqr = — 1) replaces uniform rotation: as a consequence the system is overstable. It must 
be noted that the growth rate of the amplitude of the overstable mode shown in panel 4b is not negligible: the perturbation, 
which is assumed of the order of 1% at r = 0, becomes nonlinear at r ~ 25. Model 4c also rotates differentially with Fur — — 1, 
as model 4b, but now the value of the azimuthal wave-number is smaller by a factor of 50 than in model 4b, meaning that the 
perturbation is nearly axisymmetric: this mode, which is relevant for disturbances very elongated along <j>, is clearly unstable. 

(iv) Case with Hi c > and cD t h < 0. In this case both thermal conduction and cooling are present. Let us focus again on 
cases in which 7' < 7 and TnR > —2 (outward-increasing specific entropy and angular momentum). Considering the initial 
value of the dissipative frequency cj d (0) = Q,{Hi c k' 2 + uj t h), where k' 2 = Ic'r + k! 2 + 1 = fc 2 (0), we distinguish different cases 
according to the sign of u>d(0). When cj d (0) > 0, all explored cases of non-axisymmetric perturbations in the presence of 
differential rotation turn out to be stable, so differential rotation and deviation from axisymmetry do not have destabilizing 
effects on configurations stabilized by thermal conduction. When ojd(0) < we find stability, overstability or instability 
depending on the values of the parameters. Given that in the limit of uniform rotation or axisymmetric perturbations we 
always have instability in this parameter regime (see right-hand panels of Figs. 1 and 3), it is clear that non-axisymmetry 
of the perturbations and differential rotation can have a stabilizing effect. An example is shown in Fig. 5, plotting the time 
evolution of non-axisymmetric perturbations in three configurations (5a, 5b and 5c), having the same values parameters as 
the corresponding configurations considered in Fig. 4, but uj c k' 2 =0.1, so that the initial value of the dissipative frequency is 
Wd(0) = — 0.2fi. The uniformly rotating case (model 5a) is unstable, consistent with the results of Section 4.1. Interestingly, 
the differentially rotating case with the same wave- vector (model 5b) is stable: in this case is negative at t = 0, but 
soon becomes positive because k 2 R increases as a consequence of differential rotation (making also k 2 and uj c increase). We 
recall that the analogous configuration, in the absence of conduction is overstable (model 4b in Fig. 4). Physically, the shear 
tends to distort any overdense region by making it narrow in the R direction, so the growth of the perturbation is more easily 
damped by thermal conduction. As models 4a and 4b, also models 5a and 5b are characterized by modes with k<f, of the order 
of kR and k z , so they are of interest for blob-like perturbations. If the azimuthal wave- number is small enough (i.e. the mode 
is almost axisymmetric, relevant to overdensities very elongated along <f>) the system appears unstable even in the presence of 
differential rotation (model 5c). Formally, also in this case Ud(t) becomes positive over sufficiently long times, and eventually 
the perturbation is damped. However, this must not necessarily happen in reality, because it is clear that the linear analysis 
breaks when at least one among vr, v z and p becomes much larger than unity. 

4-2.2 Baroclinic distributions 

In the previous Section we focused on barotropic distributions, i.e. configurations in which — Q(R) and p = po(po)- We 
consider here the problem of the stability against non-axisymmetric perturbations of systems with baroclinic distributions: 
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in other words, we allow for the fact that can depend on z as well as on R (implying that po is not stratified with po; see 
Section 2). In terms of the formalism here adopted, we describe solutions of equations (64-66) in cases with Fr> z =fc 0. 

As it happens for the barotropic cases, also baroclinic configurations are found to be stable, unstable or overstable 
depending on the specific choice of the parameters. In general we find that adding a relatively small vertical gradient of 
O does not change the stability properties of a system. However, we found that in some cases the fact that Q depends on 
z can have destabilizing effects. In order to isolate this phenomenon, we considered models with the same parameters as 
the barotropic models discussed above, but with Tnz ^ 0. For small values of Tsiz the baroclinic models behave like the 
corresponding barotropic models, but — in the absence of thermal conduction — instability replaces overstability for sufficiently 
large \Tq z \ (i.e. sufficiently strong dependence on z of the angular velocity). An example of thermally unstable baroclinic 
configuration is model 4d in Fig. 4, which has the same values of the parameters as the (overstable) barotropic model 4b, but 
Tnz = —0.6. When thermal conduction is present baroclinic models are found to be stable if the "corresponding" barotropic 
models are stable: compare, in Fig. 5, the barotropic model 5b with the baroclinic model 5d. 



5 IMPLICATIONS FOR GALACTIC CORONAE 

Here we use the results obtained above to address the question of whether the coronae of disc galaxies can fragment, via 
thermal instability, into cool, pressure-supported clouds, which, in the case of the Milky Way, would be identified with the high- 
velocity clouds. As already pointed out, the physical properties of the galactic coronae of disc galaxies are poorly constrained 
observationally. Such coronae are expected to be stratified, almost in equilibrium at the system's virial temperature, similar to 
the hot atmospheres of massive elliptical galaxies and galaxy clusters, but characterized by significantly lower gas temperature 
and density, and possibly by non-negligible rotation. 

MRB studied the problem of the thermal stability of the hot atmospheres of galaxy clusters, using Eulerian plane-wave 
perturbations and assuming that the gas does not rotate, and concluded that these systems are not prone to significant 
thermal instability. It has been pointed out that when a background flow is present, the use of Lagrangian perturbations is 
preferable to that of Eulerian perturbations (Balbus 1988). However, the main results of MRB have been confirmed by the 
Lagrangian study of Balbus & Soker (1989). BNF applied the thermal-stability analysis of MRB to non-rotating models of the 
coronae of disc galaxies, finding that the result obtained for the highest-temperature atmospheres of clusters (thermal stability 
or overstability) apply also to these lower-temperature systems, unless the gas distribution is perfectly adiabatic (which is 
unexpected). In the present paper we tried to verify whether the conclusions drawn from the analysis of the non-rotating 
cases extend also to rotating coronae. As already stressed, for simplicity we treated the gas as unmagnetized, so our results 
apply to real systems only to the extent that their magnetic fields do not influence substantially the behaviour of thermal 
perturbations. 

From the calculations reported in the above Sections it is apparent that rotation introduces some mathematical complexity 
in the stability analysis. As a consequence, it is hard to draw general conclusions about the thermal-stability properties of 
rotating coronae, which can be thermally stable, overstable or unstable, depending on the distribution of specific entropy 
and angular momentum, but also on the nature of the perturbations. If we knew the density, temperature and specific 
angular momentum distributions of the corona of the Milky Way (or of an external disc galaxy), the dispersion relations 
and the differential equations derived above could be used to determine unambiguously whether a given perturbation grows. 
Unfortunately, the uncertainties on the physical properties of the coronae are such that we cannot be conclusive about their 
thermal stability or instability. Nevertheless, some general considerations can be done. 

First of all, our calculations show that the stabilizing effect of thermal conduction found in non-rotating systems occurs 
also in the presence of rotation (see cases with w d > in Sections 3 and 4). As in the non-rotating case, thermal perturbations 
with sufficiently short wavelength do not grow. The interesting question is then whether perturbations with sufficiently long 
wavelength (i.e. small enough wave-number k, so that uid < 0) grow. It seems reasonable that galactic coronae are characterized 
by outward-increasing specific entropy 2 and angular momentum, so in the present discussion we can specialize to configurations 
with wfjy > and w 2 ot > 0. In Section 3 we showed that when uj d < 0, such configurations are unstable against axisymmetric 
perturbations, in contrast with analogous non-rotating configurations, which have been shown to be overstable. Taken at face 
value this result indicates that the presence of even slow rotation can modify qualitatively the stability properties of a corona: 
formally, we should conclude that rotating galactic coronae are thermally unstable, because there is at least one mode that 
grows (in fact, all the axisymmetric modes grow). However, such instability must not necessarily have important consequences 
in the practical problem of the thermal stability of galactic coronae. We note that the high-velocity clouds, which would be 
the end-products of the thermal instability of the Milky Way corona, do not appear axisymmetric with respect to the rotation 

2 BNF showed that, at least in the absence of rotation, very special conditions must be met in order to have a corona with flat specific 
entropy profile. 
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axis of the Galaxy. Moreover, axisymmetric structures of cold gas far from the plane are not observed, in general, in external 
disc galaxies. 

It is then interesting to study the behaviour of non-axisymmetric perturbations, which are relevant to the problem of 
whether blob-like, cool clouds can condense out of the corona. We have shown in Section 4 that at least some differentially 
rotating configurations, which are unstable against axisymmetric perturbations, are overstable (model 4b in Fig. 4) or stable 
(model 5b in Fig. 5) against non-axisymmetric perturbations with sufficiently high azimuthal wave-number (specifically, to 
modes with k$ of the order of ha and k z , which describe blob-like disturbances). It follows that an overdense, cool blob in 
a differentially rotating corona is unlikely to condense, even if the thermal-instability frequency is faster than conductive 
damping (i.e., cjd < 0). 

From our analysis it also emerged that the distribution of specif angular momentum is particularly important for the 
thermal stability of the coronae. For instance, uniform rotation (unexpected in real systems) would favour thermal instability, 
while (more realistic) differential rotation with outward-increasing specific angular momentum tends to stabilize. A vertical 
gradient of Q (fi decreasing for increasing z) , might be expected if the kinematic properties of the coronae reflect those of the 
neutral extra-planar gas of disc galaxies (e.g. Fraternali 2009, and references therein). In the absence of thermal conduction 
such a vertical gradient, if sufficiently strong, might have a destabilizing effect (compare baroclinic and barotropic models 
in Fig. 4). However, such an instability is not expected to occur in the presence of even highly suppressed conductivity. In 
fact, the plots in the right-hand panel of Fig. 5 show clearly that conductive damping is enhanced by the shear induced by 
differential rotation: even non-axisymmetric perturbations that have initially long enough wavelength — such that Wd(0) < — 
are effectively sheared and therefore stabilized by thermal conduction, independently of whether the distribution is barotropic 
(model 5b) or baroclinic (model 5d). 

Altogether these considerations lead to the conclusion that, generally speaking, differential rotation does not make galactic 
coronae thermally unstable to non-axisymmetric perturbations. We cannot exclude that thermal perturbations grow in specific 
cases, but a definitive answer to the question of whether the corona of a given galaxy is thermally unstable can be obtained 
only when a detailed model of the corona, including its specific angular momentum distribution, is available. 



6 SUMMARY AND CONCLUSIONS 

Motivated by the question of whether cool, pressure-supported clouds can condense out of the hot coronae of disc galaxies, 
we studied the problem of the thermal stability of rotating stratified fluids, in the presence of radiative cooling and thermal 
conduction. As in the non-rotating case, the time evolution of a perturbation depends on its wavelength, so it is useful to 
distinguish short-wavelength perturbations (such that the dissipative frequency cj d is positive) from long-wavelength pertur- 
bations (such that the dissipative frequency ujd is negative) . We found that — against either axisymmetric or non-axisymmetric 
perturbations — a uniformly rotating, convectively stable configuration is thermally stable when w d > (damping by thermal 
conduction), but is thermally unstable when aid < 0. Similarly — against axisymmetric perturbations — a differentially rotat- 
ing, convectively stable corona with outward-increasing specific angular momentum is thermally stable when tOd > 0, but is 
thermally unstable when cjd < 0. Non-axisymmetric perturbations in the presence of differential rotation behave differently 
from the axisymmetric ones. In the absence of thermal conduction, the combination of non-axisymmetry of the perturbation 
and of differential rotation has a stabilizing effect (turning instability into overstability) , and barotropic distributions tend 
to be more stable than baroclinic distributions. In the presence of thermal conduction, stability replaces overstability: in 
differentially rotating systems, the shear makes conductive damping of non-axisymmetric disturbances particularly effective 
for both barotropic and baroclinic distributions. 

These results have been discussed in the context of the problem of the growth of thermal perturbations in galactic coronae. 
The above calculations allow to address unambiguously this question when applied to a specific model of rotating galactic 
corona. Unfortunately, given that the physical properties of the hot atmospheres of disc galaxies (in particular, the distribution 
of specific entropy and angular momentum) are very poorly constrained observationally, it is difficult to be conclusive about 
whether these systems are prone to thermal stability. Additional uncertainties come from magnetic fields, which are expected 
to be present in real systems, but have been neglected for simplicity in the studied models. The question of the effect of 
magnetic fields on the thermal instability is very interesting and would require a full thermal-stability analysis of a rotating, 
magnetized corona, for which our calculations may be a starting point. 

Though limited by the mentioned uncertainties, the present work gives further indications against the hypothesis that 
thermal instability is important for galactic coronae. In particular, though rotation can destabilize against thermal disturbances 
very elongated in the azimuthal direction, we argued that blob-like thermal perturbations are unlikely to grow in a differentially 
rotating corona. This finding, combined with previous results on non-rotating coronae (BNF) , suggests that the high- velocity 
clouds of the Milky Way did not form spontaneously from small thermal perturbations in the Galactic corona, but must be 
of external origin. A possibility is that at least the seeds of these clouds are formed originally as a consequence of stripping 
of gas-rich satellites or cosmic infall of cold gas. These cool gaseous seeds might then grow while travelling through the hot 
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Galactic corona (Sommer-Larsen 2006; Keres & Hernquist 2009), for instance via turbulent mixing (Marinacci et al. 2010). 
Provided that the overdensities associated with these accreted seeds are substantial (i.e. they are non-linear perturbations), 
such a scenario for the formation of the high- velocity clouds is not in contrast with the results of the present paper, in which 
only linear perturbations are considered. 
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APPENDIX A: SHORT- WAVELENGTH, LOW-FREQUENCY, AXISYMMETRIC PERTURBATIONS 
ARE ALMOST ISOBARIC 

Here we show that in a differentially rotating fluid short-wavelength, low-frequency, axisymmetric perturbations are almost 
isobaric, in the sense that p/po <C p/po, where po and po are the unperturbed pressure and density, while p and p are the 
pressure and density perturbations. Let us consider the linearized mass and momentum equations (9-12): eliminating vr, v z 
and v<f>, and we obtain 



(Al) 



p a p 
po P po ' 
where 

il) 2 .u) 2 k R A pR £IQ Z knAp Z .k z A pz 
a = W~ l ^^^~ ~ (A2) 
and 

a _ ^ 2 fcfl , k R k z k'j 

P = ^^ +2 ~~k 2 ~ + k^' [AS) 
with uj 2 = Co 2 — 2Q£Lr — 2Q, 2 . In the limit of short wavelengths and low frequencies adopted in Section 3, all terms of a are 
infinitesimal, while f3 is finite. It follows p/po <C p/po and then the perturbed equation of state p/po = p/po — T /To can be 
approximated as p/po ~ — T /To (see Tribble 1989, for the analogous calculation in the non-rotating case). 



APPENDIX B: SOME PROPERTIES OF CUBIC EQUATIONS 

We recall here a few properties of cubic equations. Let us consider the equation 
h 3 + an 2 + bn + c = 0, 

in the unknown h, with coefficients a, b and c. The associated discriminant is 

A = -4a 3 c + aV - 46 3 + 18abc - 27c 2 . 

The roots of equation (Bl) have the following properties: 



(Bl) 
(B2) 



- Number of real roots. If A > the equation has three real roots; if A < the equation has one real roots and a pair of 
complex conjugate roots. 

- Routh-Hurwitz Theorem. All the roots have negative real parts if and only if the following conditions are satisfied: 



a > 0, 



a 1 

c b 



>0, 



a 1 
cb a 
00 c 



> 0. 



- Viete's formulae. The three roots hi, fi2 and fi3 satisfy Viete's formulae: 

hi + hi + hz = —a, 
hihi + hihz + hzhi = b, 
hihihz = —c. 
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